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The present uncertainties in the knowledge of the neutrino cross sections for E v ~ 1 GeV, that is in the energy 
range most important for atmospheric and long baseline accelerator neutrinos, are large. These uncertainties do 
not play a significant role in the interpretation of existing data, however they could become a limiting factor in 
future studies that aim at a complete and accurate determination of the neutrino oscillation parameters. New 
data and theoretical understanding on nuclear effects and on the electromagnetic structure functions at low Q 2 
and in the resonance production region are available, and can be valuable in reducing the present systematic 
uncertainties. The collaboration of physicists working in different subfields will be important to obtain the most 
from this available information. It is now also possible, with the facilities developed for long baseline beams, to 
produce high intensity and well controlled i/-beams to measure the neutrino interaction properties with much 
better precision that what was done in the past. Several projects and ideas to fully exploit these possibilities are 
under active investigation. These topics have been the object of the first ^-interaction (NUINT01) workshop. 



1. Introduction 

Neutrino physics is living in a "golden era" . 
Experiments with atmospheric and solar neutri- 
nos have given strong evidence for the existence of 
flavor oscillations, or (less likely) some other form 
of 'New Physics' beyond the Standard Model 
jl). The study of the v flavor transitions offer 
the possibility to obtain information about the 
neutrino masses and mixing, and the knowledge 
about these quantities hopefully represents a fun- 
damental window on the physics of the unifica- 
tion. Theorists are busy trying to make sense of 
the unexpected results that have been obtained 
so far, on the other hand new experimental stud- 
ies are planned to measure with greater precision 
the properties of the v flavor transitions. These 
future studies, are designed to determine if the v 
flavor evolution is completely described by the 3- 
flavor mixing model, or if a more complex dynam- 
ics (transitions to additional sterile states, neu- 
trino decay, flavor changing neutral currents) is 
necessary, and should measure the entire set of 
parameters (2 squared mass differences, 3 mixing 
angles and one CP-violating phase) that deter- 
mine the v flavor evolution in space-time in the 
3-flavor scenario. 

The discovery of v oscillations has been per- 



formed with natural (solar and atmospheric) neu- 
trino sources, however most of the future stud- 
ies will be performed with artificial neutrinos. In 
particular, long baseline (LBL) accelerator neu- 
trino beams will play a fundamental role. The 
first of these projects, the KEK to Kamioka 
(K2K) z/-beam, with (E u ) ~ 1.8 GeV and L ~ 
250 Km, has already collected data, that support 
the existence of the — > v T transitions as in- 
dicated by the atmospheric neutrino data 
Other projects are already funded and in con- 
struction (Fermilab to Minos ||[| , and CERN to 
Gran Sasso |5U|], both with L ~ 730 Km), or in 
the design phase (the JHF project |0,||). Sev- 
eral 'next generation' projects are also in discus- 
sion, involving "super-beams" , (very intense v— 
beams produced with the traditional technique: 
p + Target — > 7r ± , K^, . . . — » v^{p^)), or neutrino 
factories (muon storage rings that produce ex- 
traordinarily well controlled neutrino beams from 
muon decay: /i + — > 77 M + v e + e + ). 

The v energy (and pathlcngth) in the LBL 
projects is determined by their scientific goals 
and the existing technical and experimental con- 
straints. The range of energy E v ~ (0.2-10) GeV 
will have crucial importance. In fact the differ- 
ence in quantum mechanical phase developed by 
two mass eigenstates after propagation for a dis- 
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Figure 1. Top panel: Energy distribution of the 
event rates (dN^/dln E v ) for atmospheric neutrinos 
and for operating (K2K) or planned long baseline 
neutrino experiments [Fermilab to Minos (LE) and 
CERN to Gran Sasso]. Bottom panel: charged cur- 
rent Vfj, cross section calculation || compared with 
data. Note the importance of quasi-elastic and single 
pion production. 



tancc L is: 

A<p jk = (E k - E s ) L ~ (Am% L) / (2E U ) (1) 

where Am 2 fe = m 2 — m 2 is the v squared mass 
difference. The probability for flavor transitions 
becomes important when at least one quantum 
phase difference is of order unity, and this de- 
termines the combinations of L and E v that are 
needed to study the flavor transitions once the 
order of magnitude of the Am 2 are known. The 
largest squared mass difference, measured in at- 
mospheric v experiment: |Am 2 tm | = jAm 2 ^! is 
of order 3 x 10 -3 eV 2 . A quantum mechanical 
phase difference Atp 23 = tt corresponds to the en- 
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At this energy one expects a maximum of the 
Ufj, — » v r oscillation probability for a distance L 
between v source and detector. Note that also the 
evidence for oscillation for atmospheric neutrinos 
has been obtained with neutrinos of energy E v ~ 
1 GeV. The top panel of fig. 1 shows the energy 
spectra of different long baseline projects, and for 
atmospheric neutrinos^]. 

The flavor composition, and energy spectrum 
of a v beam in a detector (located at a distance 
L from the source) is determined from the obser- 
vation of neutrino interactions, and clearly un- 
certainties in the knowledge of the neutrino cross 
sections are a limiting factors in the sensitivity. 
There are uncertainties about the absolute value 
of the cross sections, the ratio cr nc /a cc between 
the neutral and charged current interaction rates, 
the average multiplicity in the final state, the en- 
ergy and angular distribution of final state lepton 
(and hadronic particles) . Depending of the detec- 
tor type and the scientific question asked, these 
uncertainties can have a different importance. 

The NUINT01 workshop^ has been dedicated 
to a discussion of the problem of the neutrino in- 
teraction properties, attempting to estimate: (i) 
the size of the existing uncertainties, (ii) their 
probable impact on studies of the fundamental 
v properties, and (iii) possible methods to reduce 
these uncertainties. A list of relevant questions 
is: 

(A) How well have we measured the v cross sec- 
tions and interaction properties? 

(B) How good is our theoretical understanding 
(and computational capabilities) of the v 
cross sections and interaction properties? 

(C) How well do we need to know u v to deter- 
mine with the desired accuracy the v flavor 
evolution and mass matrix? 



1 The CERN Gran Sasso project is designed to study the 
appearance of r neutrinos, with the process — > v T — * 
t~ . Since the threshold for the r production is ~ 3.5 GeV, 
the beam is designed to have a larger average v energy. 
2 See the web page neutrino.kek.jp/nuint01/ 
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(D) What experimental programs can/should 
be planned in order to improve our knowl- 
edge and understanding of the v interaction 
properties? 

2. Measurements of ct„ 

The neutrino cross section can be obtained 
from data, from the formula: 



<J V {E V ) 



(3) 



where A^ nt is the number of detected interactions, 
<&„ is the neutrino flux, and £dct is detector effi- 
ciency. This measurement has been performed in 
the energy region of interest [E v ~ 1 Gev) by 
several experiments Jl^]-Jl^]. The most accurate 
one were performed with bubble chambers filled 
with hydrogen or deuterium. A review of this 
data is given in fIo|| . The uncertainty of <j v is 
dominated by systematic errors in the determi- 
nation of the v flux. For example for the data 
of [[llj (the 12-foot hydrogen/deuterium bubble 
chamber at the Argonne National Laboratory), 
the authors estimated an overall 15% uncertainty 
for the neutrino beam that was produced by a 
12.4 GeV proton beam interacting on a Beryllium 
target. The systematic uncertainty was estimated 
as as a convolution of different factors: uncer- 
tainties on the fundamental 7r production cross 
section (5%), uncertainties of pion reabsorption 
in the target (5%), proton beam intensity (2%), 
and smaller contributions. Recent critical analy- 
sis jl7| of the uncertainties in the p-nucleus cross 
sections, related to the calculations of the atmo- 
spheric neutrino fluxes have concluded that the 
uncertainties on the hadronic cross sections are 
larger (of order 15-20%). To be conservative, it 
could be safer to consider a systematic larger that 
what was quoted in the original experiments. The 
facilities developed for the long-baseline experi- 
ments can provide v beams significantly more in- 
tense that the beams used in the 70 's. In order 
to use these beams to obtain absolute values for 
the cross sections it is essential to have good con- 
trol of the normalization and energy shape of the 
neutrino beam. 



The K2K experiment has already collected in 
the near detector a very important sample of 
events (see section 4), and a problem of consid- 
erable interest is to determine how well the v 
beam is known. The same question clearly applies 
for the near detector in the Fermilab to MINOS 
project. 

In principle it is possible to determine the 
quasi-elastic cross section on (quasi)-free nucle- 
ons without a precise knowledge of the v flux, 
making use of some theoretical input. If the in- 
cident v energy is known, from a measurement 
of the energy and direction of both final state 
particles in the scattering v^n — > pTp one can 
determine both the neutrino energy E v and the 
Q 2 of the reaction^ It is therefore possible to de- 
termine the shape of the differential cross sec- 
tion da I dQ 2 (Q 2 ) , that to a good approxima- 
tion can be considered as a function of a sin- 
gle parameter the axial mass Ma (see next sec- 
tion), while the absolute normalization, that is 
the value of the cross section at Q 2 — ► can 
be inferred from (3 decay measurements. This 
program has been performed in [ fl2[ (data from 
the 7-foot hydrogen/deuterium bubble chamber 
at the Brookaven National Laboratory) with the 
result M A = 1.07 ± 0.06 GeV, and in @ (ANL 
12— ft bubble chamber) with the result Ma 
1.01 ± 0.09 GeV. The total cross section (and 
therefore the event rate) is to a good approxi- 
mation proportional to M\. In |ll| the analy- 
sis of the event rate results in a lower estimate 
Ma = 0.75 ± 0.12. The poor agreement between 
the two methods suggests the existence of system- 
atic effects. 

Several other experiments have measured v in- 
teractions in the energy range of interest using 
nuclear targets. In the published analysis of these 
works very often the results are corrected for nu- 
clear effects, and the data are the expressed as 
cross sections on free nucleons. Since the nu- 
clear effects themselves are difficult to compute 
(as discussed in the following) and the subject of 
some controversy, some careful considerations are 
needed to use these results. For practical pur- 



3 The v reaction has to be performed on a nuclear target 
(preferably deuterium) and nuclear corrections have to be 
considered. 
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poses it clearly would be very valuable to have 
available also the results for interactions on nu- 
clear targets, before the nuclear corrections. 

An example of this problem can be seen 
analysing the results of the two v experiments Jll| 
[ig[| operating with the v-beam at the Institute for 
High Energy Physics (IHEP) at Serpukhov (Rus- 
sia), that were presented at this workshop. The 
SKAT bubble chamber jl8| has collected data on 
a set of low multiplicity reactions, namely the 
final states (p,~p), (/i~p7r + ), (/i~p7r°), (/i~n7r + ) 
with an additional non-observed 'nuclear specta- 
tors' state. Each one of the considered final states 
can be obtained as the result of nuclear cascading 
processes, such as ir absorption (tt + A — > A*), tt 
production (pp — > pnir + , pp — » ppw° , . . .); and tt 
charge exchange (for example tt + ti «-> n°n). The 
experimentalists have calculated the probabilities 
for all these processes with montecarlo methods 
to estimate the rate for exclusive reactions on free 
nucleons, however these corrections are probably 
the largest uncertainty in the problem. 




W (GeV) 



Figure 2. Cross section for the scattering + n — » 
H~ + X. The thick red line is a calculation that uses 
the DIS formula with the structure functions given 
by the naive quark-parton model and the PDF's of 
GRV-94 pi] ]. The other lines are the cross sections 
for the exclusive production of a resonance (y^ + n — > 
fj,~ + R) and their sum, according to the Rein and 
Sehgal model E§. 



3. Cross section Formulae 

Only the interactions of neutrinos with nucle- 
ons and nuclei are relevant for oscillation studies, 
since the ratio a ue / '(TupOn) 1S °f order m e /m p ; the 
cross section can naturally be separated into a 
charged current (CC) and a neutral current (NC) 
component: <J Va = cr£ c + er"^. For low energy 
neutrinos it is also natural to decompose the cross 
section according to the multiplicity: 



a Ua (qel) 



;(1tt) +o-(2tt) 



The correct calculation of the quasi-elastic and 
one pion components has great importance in the 
analysis of existing data since most of the data 
on atmospheric neutrinos has been obtained with 
water Cerenkov detectors (Super-Kamiokande, 
and before 1MB and Kamiokande) that, because 
of relatively poor pattern recognition capabilities, 
have limited most of their analysis to the so called 
single-ring events, where only one particle is vis- 
ible. This selects quasi-elastic scattering^ with 
in addition events with one (or more) charged 
pion(s) below the Cerenkov threshold. The bot- 
tom panel of fig. 1 shows a plot of separated 
according to the pion multiplicity and compared 
with some existing data. The importance of the 
low multiplicity channels is clearly displayed 

3.1. Quasi elastic scattering 

The quasi-clastic cross section can be written 
in terms of nucleon form factors. The most gen- 
eral matrix element for a nucleon transition can 
be written j20| in terms of six form factors: 

{p{p')\J a {Q)\n{p)) =u p (p') 



F v ( g 2 ) 7 « + zF M ( g >^ 
F A (q 2 h a 7 5 + F P (q 2 ) 2Af 



F v s ' 



75 



(l 2 ) q TT 



in(p) 



The hypothesis of CVC allows to relate the vec- 
tor form factors Fy and Fm to the electromag- 
netic form factors of the proton and neutron; the 
contribution of second class currents (that violate 
G-parity) can be assumed to vanish, and one re- 
mains with two unknown functions the axial form 



4 Most final state protons are below the Cerenkov thresh- 
old. 
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factor Fa and the pseudo-scalar form factor Fp. 
The contribution of the pseudoscalar form factor 
to the cross section is proportional to m 2 /M 2 , 
therefore this term is important only for v T , while 
the axial form factor Fa remains undetermined. 
The value of Fa at Q 2 = can be related to the 
axial coupling measured in /3-decay experiments 
(F A (0) = 9A = 1-2670 ± 0.0035), but the Q 2 de- 
pendence of Fa has to be determined experimen- 
tally. An often used parametrization is a dipolar 
form: 



Fa(Q 2 ) = 9A 



1 



(4) 



that depends on the single parameter Ma- There- 
fore in first approximation the problem of the de- 
termination of the quasi-elastic cross section can 
be indcntificd with the measurement for MA. A 
review of the status of our knowledge of the nu- 
cleon form factors was given in |35[ and p6| . 

3.2. Inelastic channels 

The kinematics of a final state lepton is de- 
scribed by two variables, that correspond to its 
energy Eg and angle with respect to the initial 
neutrino 0£ U . The kinematical variables can be 
chosen as relativistic invariants such as the adi- 
mensional scaling variables x and y; of particular 
dynamical importance are the transfer momen- 
tum in the interaction: quantities: 



Q 



-(f 



\Vv -pif 



(5) 



and the square of the hadronic mass in the final 
state: 



W 2 =p 2 ad = ( Pl + q) 2 



(0) 



The kinematically allowed range of Q 2 is the in- 
terval between Q^ in — and Q, 2 nax — 2m p E„ 
that grows linearly with energy. For Q, 2 nax <C 
M 2 , (that is for E v <C 3 TeV) the effects of 
the W propagator are small and the cross sec- 
tion a v grows approximately linearly with en- 
ergy. Similarly the allowed range of the squared 
hadronic mass W 2 is between the values m 2 and 



,E V . 



2m, 

In general [|37jJ3q] , the inclusive differential neu- 
trino cross section can be written in terms of 5 



structure functions Fj(x,Q 2 ), according to the 
general formula: 



dxdy 
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where the two signs for the F3 component cor- 
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Figure 3. Calculated distributions of the differen- 
tial cc event rates {dN^/ dW)'m the far detector as a 
function of W the hadronic mass in the final state. 
The top panel is for the K2K experiment, the bottom 
panel for the Fermilab to Minos project. The vertical 
line at W — Mn represents quasi elastic scattering, 
the dashed line is calculated according to the DIS 
formula, the solid line is the contribution of the reso- 
nances according to the Rein and Sehgal model [E3|. 
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respond to v and V. The scaling variables x 
and y are denned as y — (q ■ Pi)/{p v • Pi) and 
x = Q 2 /(2iWjv E v y). Note that the structure 
functions F4 and contribute terms propor- 
tional to m^/(2Mjv E v ) and are therefore of prac- 
tical interest only in the case of v T cc interactions 
(where they constitute a source of uncertainty). 
It is well known that in the Deep Inelastic Scat- 
tering (DIS) limit, that corresponds to Q 2 and 
W 2 large, the structure functions can be related 
to the Parton Distribution Functions that can be 
well measured in different experiments, and have 
a calculable Q 2 evolution. For example in leading 
order of an a s expansion: 

F%' cc = 2xF^ cc = 2x(d + s + u + c) 

F^ cc = 2(d + s~u-c) ( '> 

Therefore in the DIS region the neutrino cross 
section is well understood and under good the- 
oretical control. However, essentially the entire 
phase space relevant for atmospheric neutrino 
contained events, is well outside of the DIS re- 
gion, and this remains true for the K2K, Minos 
to Fermilab and JHF projects. An illustration of 
this problem is shown in fig. 2, where we plot the 
cross section for the scattering a u +n — > +X &t 
an energy of 3 GeV. At this energy (that is actu- 
ally rather high for atmospheric v and most LBL 
experiments), the range of possible W is small 
and extends in the region where hadronic reso- 
nances are present. The naive use of the DIS 
formula with the PDF's of Gluck Reya and Vogt 
]2l[ | results in a smooth curve, that predicts a 
non vanishing cross section also in the interval 
m p < W < mp + m^. Obviously a realistic calcu- 
lation must predict a more complex structure, as 
it is the case for the model of Rein and Sehgal |2lj| ] 
whose results are also shown in the figure. Fig- 
ure 3 is another illustration of the importance of 
the resonance region in the K2K and Fermilab to 
MINOS projects. The delta function for W — m p 
has an area of ~ 42% for K2K and of order 11% 
for MINOS (LE beam). 

A central problem for the calculation of the v 
cross section is the correct description of the low 
Q 2 , low W 2 kinematical range, where perturba- 
tive QCD is not valid. Luckily, in recent times, 
a large body of highly accurate data on electron- 
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Figure 4. Structure function F| m ' as a func- 
tion of x for four values of Q 2 , with a logarithmic x 
scale. The symbols indicate various experiments. The 
solid curves in the top two panels represent the calcu- 
lated distributions from the GRV collaboration pi) , 
evolved from Q 2 = 0.4 (GeV/c) 2 . The solid curve 
in the third panel represents the input distribution 
itself. From @. 



proton scattering in this kinematical region has 
become available, mostly from the Jefferson lab- 
oratory (see |39|- |4l|j38[ ) . An example of this data 
is shown in fig. ^ where is plotted the electro- 
magnetic structure function F^ ™ 1 that in the DIS 
limit in leading order can be expressed in terms 
of PDF's as: 

FS- m -=T,e 2 j x[q J {x)+^ j ] (8) 
j 

One can see that the production of resonances 
plays a very important role, and that three mass 
peaks are clearly visible at low Q 2 . Fig. || also 
illustrates the importance of a proper treatment 
of the relevant effects in the low Q 2 , region. The 
physical insight that can be obtained from the 



7 



0.0 





1 1 1 1 1 


1 ' ' ' ' 1 ' ' ' ' 1 ' ' ' ' - 
Quark- partem model: PDF = GRV-94 - 






Q 3 = 30, 10, 3, 1, 0.4 GeV £ — 






V \0.4 — 




30\ 






1 , , , , 





0.0 0.2 0.4 0.6 0.8 1.0 

X 



Fits to e-p data 




0.0 0.2 4 0.6 0.8 1 

X 



Figure 5. Structure function F| m (a;,Q 2 ) for e- 
p scattering. The top panel shows what is obtained 
from the simple use of the naive quark-parton model. 
The bottom panel shows fits to the SLAC and JLab 
data. 



electron scattering data can be applied also to 
neutrino interactions, even if significant theoreti- 
cal uncertainties remain. A concept of great im- 
portance in the resonance region is the notion 
of "parton-hadron" duality. Already in 1971, 
Bloom and Gilman |2?]] observed empirically, that 
the scaling structure function F% m (x) that char- 
acterizes electromagnetic scattering in the DIS 
(high Q 2 region), represents a good average over 
the resonance bumbs of the inclusive data ob- 
served at much lower Q 2 . The agreement can 
be significantly improved rescaling the x vari- 
able for target mass correction and using for 
example the Nachtmann variable £ = 2a;/ (1 + 
v /l + 4M 2 x 2 /Q 2 ). 

3.3. Nuclear effects 

Neutrino detectors are usually composed of 
heavy nuclei (the most important nuclear species 



are oxygen (SK), argon (Icarus), iron (Minos) 
and lead (Opera)), and in the energy range we 
are discussing nuclear effects are important and 
have to be taken properly into account. Schemat- 
ically one can separate the nuclear effects into 
two parts: (i) the v interaction with a bound nu- 
cleon, and (ii) the propagation of the final state 
hadrons in the nuclear matter, even if a clear dis- 
tinction of these two effects is not rigorously valid. 
The simplest and most commonly used approxi- 
mation to describe the interaction with a bound 
nucleon is the Fermi gas model. In this model 
the nucleus is described as an ensemble of nucle- 
ons with a momentum distribution f Pt n(Pi)- The 
energy that corresponbds to the 3-momentum pi 
is Ei — \Jpf + M 2 — Eb where the quantity £s is 
the binding energy that can be different for pro- 
tons and neutrons^. In the simplest version of the 
model p^ j, the initial nucleon momentum distri- 
bution has a "zero temperature" shape: 

/(&) = e( PF - \pi\) (9) 

that vanishes when \pi\ is larger than the Fermi 
momentum pp, in other models it has a more 
complex shape that includes higher momentum 
tails. Quasi-elastic scattering can be simply de- 
scribed as an elementary (2 — > 2) process such 
as i/£ + n* — » t~ + p* where n* and p* are off- 
shell initial final states nucleons with 3-momenta 
Pi and pi + q (q — p v — pt is the transfer 4- 
momentum) and binding energies Ei and e f . This 
procedure allows naturally to implement Pauli- 
blocking, forbidding the scattering into occupied 
nucleon states. This can be done including a fac- 
tor [1 — f p (pi + q)] in the calculation of the matrix 
element for the scattering. Some natural ques- 
tions, that have been discussed at the workshop 
are: (i) how good is the Fermi gas model as a de- 
scription of the nuclear effects; (ii) how large are 
the expected corrections to the Fermi gas model; 
(iii) what is the best (or most convenient) theoret- 
ical framework to describe nuclear effects beyond 
the Fermi gas model. It is clear that conceptually 
the Fermi-gas model is a very simple and non- 
realistic description of the scattering |24|,|2^] , it is 

5 In some implementation of the model eb is also a function 
of pt. 



8 



,6 (v,n )'V 




0.10 0.20 0.30 0.40 



•a 




0.40 0.60 
E^IGeV] 



0.80 1.00 



Figure 6. Cross section integrated on the muon 
emission angle as a function of the muon energy. The 
full lines have been calculated with the Fermi gas 
model (FG label) and the dashed lines with the model 
described in that includes final state interactions 
(FS1 label), (from @). 



however not entirely clear how large is the error 
introduced by this simplified treatment. Before, 
during and after the workshop, it was attempted 
to produce and compare different "beyond the 
Fermi-gas model" calculations, (see the critical 
discussion by R. Seki E| in these Proceedings). 
Hopefully a generally accepted model will soon 
emerge from this work. An example of the re- 
sults for one of these more sophisticated models 
(by G. Co and collaborators g§). is shown in 
fig. ^| This model predicts a suppression of the 
cross section by approximately 15% with respect 
to the Fermi gas model, and a softer spectrum for 
the final state lepton. 

The uncertainties due to nuclear effects become 
larger when discussing the resonance production 
processes and the propagation of hadronic parti- 
cles in the nuclear matter. It is important to no- 



tice that because of pion reabsorption processes 
the calculation of pion-less final states must also 
include contributions from resonance production. 

4. The K2K near detector data 

The K2K experiment has collected data with 
the neutrino beam obtained from a 12 GeV pro- 
tons beam at the KEK accelerator. The goal of 
the experiment is to study the existence of neu- 
trino flavor transitions comparing the event rates 
at a near detector, located approximately 300 m 
from the pion production target, with event rate 
at a far detector(Super-Kamiokande) located at 
a distance of 250 km. The event rate recorded at 
the far detector has been measured to be smaller 
that the rate extrapolated from the near detec- 
tor results, with a suppression that is consistent 
with the <-> v T oscillation interpretation of the 
atmospheric neutrino data [0J|] . 

The neutrino interactions collected at the near 
detector represent a very large sample with poten- 
tially extremely useful information. The discus- 
sion of this data has been perhaps the highlight 
of the workshop. 

The near detector for the K2K experiment is re- 
ally a combination of three detectors. The main 
flux measurement is based on a one-kiloton water 
Cerenkov detector. This detector uses the same 
photo-multipliers as Super-Kamiokande, and has 
essentially the same design (except for scale) and 
the analysis for the two detectors can use the 
same reconstruction algorithms. A fine grained 
detector sits downstream of the water detector. 
It is made as a scintillating fiber tracker with 
water targets enclosed between layers of tracking 
material. The next component of the near detec- 
tor is a stack of iron plates interleaved with drift 
tubes to serve as a muon range detector. Both 
the fine-grained and muon range detectors have 
transverse pairs of tracking planes so tracks can 
be reconstructed in three dimensions. 

The data of the three near detectors has been 
discussed in a set of 4 contributions by Ishida, 
Walter, Mauger and Mine @ ||]. The event 
sample collected up to now in the near detec- 
tors amount to a very significat statistics: 27,000 
events in the water detector, 8,300 in the fine 
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grained detector and 125,000 in the iron detector. 
A comparison of the data with a Montecarlo sim- 
ulation shows in general good agreement, demon- 
strating that there are no large systematic errors 
in the simulation chain, however small but signif- 
icant effects are visible. 

The fine grained detector |i(J as a whole is 
excellent for event classification studies since it 
can distinguish quasi-elastic and inelastic events. 
This is performed measuring the momenta and di- 
rections of the muon and the proton in events of 
type: + A — > pT + p(+ . . .). In case of quasi- 
elastic scattering on a free nucleon, because of 
4-momentum conservation, the information ob- 
tained from the the measurements is redundant, 
and the measurements of E^, 9^ and E p allow to 
determine the initial neutrino energy and the final 
state proton direction. For scattering on a bound 
nucleus, since the spectator nucleons absorb some 
4-momentum, the predicted values represent the 
central values of narrow distributions. The con- 
tribution of the inelastic channels can be sep- 
arated from the quasi-elastic one because it is 
much broader as illustrated in fig. [?]. 

A method to discriminate between the hypoth- 
esis of — > v T and — > v s t eI transitions is 
the study of the neutral current event rate. This 
event rate is suppressed in the case of transi- 
tions of into sterile neutrinos, while it re- 
mains constant for standard — » v T transi- 
tions. The most important channel for the study 
of the neutral current event rate is the reaction: 
v x + A — > v x + A + 7r° that has the lowest mul- 
tiplicity and a clear experimental signature. The 
problem is that in order to predict the event rate 
for this process one needs to know with suffi- 
cient precision the absolute value of the cross 
section, or at least the ratio with respect to a 
measurable charged current process. The K2K 
near detector is the ideal instrument to measure 
this cross section (or a cross section ratio such as 
c nc (l7'" )/o'qei- The status of these measurements 
is described in jiTj . Approximately 5000 ir° were 
reconstructed in a 50 ton fiducial volume in the 
water Cerenkov near detector, determining the 
ratio (7r°/FC M ) with small statistical error and 
an overall systematic error of order ~ 9%. This 
can be applied to the 414 reconstructed ir° in the 




Figure 7. Results from 2-track events in the K2K 
fine-grained near detector. The abscissa is the cosine 
of the angle between the second track reconstructed 
direction and the direction calculated for a proton 
based on the assumption that the other track was a 
muon and 2-body kinematics as is shown in the inset 
cartoon. The fit of this distribution can give a solid 
measurement of the non-quasi-elastic contribution in 
the data. 



atmospheric neutrino sample (1289 days) of SK 
to obtain a result that disfavors the oscillation 
into sterile neutrinos. 

It is well known that atmospheric neutrinos 
were initially considered mostly as an unpleasant 
background for proton decay studies, because for 
any detector resolution, there is a finite proba- 
bility that a v interaction will mimic a proton 
decay event. It is clearly very important to esti- 
mate correctly this probability for the interpreta- 
tion of the existing data, and more important for 
the planning of future experiments. This proba- 
bility is clearly strongly dependent on the v inter- 
action properties. The exposition of the 1 kton 
(50 tons fiducial mass) water Cerenkov detector 
to the v fluence obtained from 3 x 10 19 protons is 
roughly equivalent to a ~ 1 Mton yr exposition 
to the atmospheric z/-flux and therefore allows to 
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study experimentally the importance of the atmo- 
spheric background for p decay studies (see Q 
for more discussion). 

5. How well do we need to know the v cross 
sections? 

It is necessary to estimate carefully the possible 
impact of systematic uncertainties in our knowl- 
edge of the neutrino cross sections for the inter- 
pretation of present and future neutrino data. 

The effect of uncertatainties on a v for the in- 
terpretation of the SK atmospheric neutrino data 
has been addressed at this workshop in the con- 
tribution of K. Kaneyuki p9[ |. This work anal- 
ysed the effects of three modifications of the neu- 
trino cross section: (a) a change in ctnc of ±30%, 
(b) the setting to zero of the cross section for 
7T° coherent production on a nucleus (the process 
v+A -> v+A+tt [||,|9)), and (c) a modification 
of ±30% for the ratio CNQE/oqci- The SK simula- 
tion was repeated including these modifications, 
and the allowed region in the plane (Am 2 , sin 2 29) 
was recalculated following the same standard pro- 
cedure. The allowed intervals calculated with the 
modified cross section are remarkably similar to 
the normal ones. The effect of the ±30% increase 
of the neutral current cross section corresponds to 
a small (~ ±4%) increase in the allowed range for 
Ato 2 |. This can be understood noting that the 
neutral current events contribution to the single- 
ring events is larger for /i-like events, therefore 
a larger neutral current cross section corresponds 
to a larger (/i/e)MC, to a smaller double ratio 
(/j,/e)data/(At/e)MC an d therefore to a larger sup- 
pression of the /x-rate, and corresondingly to a 
larger |Am 2 |. 

A ±30% modification of the multi-pion cross 
section correspond to a shift of order ±0.02 in 
the minimum allowed value of the mixing param- 
eter sin 2 29. This can be qualitatively understood 
observing that in the presence of a larger contri- 
bution of the non-quasi elastic component, the 
correlation between the muon and the neutrino 
becomes weaker. The Up/Down asymmetry ob- 
served at the muon level is in general smaller than 
the true asymmetry at the neutrino level, because 
of finite angular resolution effects, with the dom- 



inant effect being the angle 9^ v between the ob- 
served muon and the parent neutrino that is de- 
termined by the interaction cross section. Assum- 
ing a weaker correlation, the inferred true asym- 
metry becomes stronger, and the mixing param- 
eter larger. 

The effect of neglecting the coherent ir° cross 
section is found to ne negligible. These results are 
very encouraging, in the sense that they demon- 
strate that the atmospheric neutrino result is very 
robust, and that the uncertainties in our knowl- 
edge of the neutrino cross section is actually very 
small. This is after all not so surprising and is a 
consequence of the fact that essentially all effects 
in atmospheric neutrinos are visible in ratios such 
as e//j, or up/down, and the theoretical uncertain- 
ties tend to cancel. 

A problem that has attracted a fair amount of 
attention in the past is the question of the ab- 
solute value of the atmospheric neutrino fluxes, 
and the precise shape of their energy dependence. 
Different calculations of the atmospheric v fluxes 
pb]-|34]| yield results that differ in absolute nor- 
malization by 20% or more, and also have non 
negligible differences in the slope of the energy 
spectrum. These differences have negligible ef- 
fects on the neutrino oscillation interpretation of 
the data, however it would be very desirable to 
determine the correct normalization of the fluxes, 
but this is of course only possible is the absolute 
value of <j v is known with sufficient precision. 

Future experiments with accelerator beams 
aiming at more precise determination of the neu- 
trino oscillation parameters will very likely have 
more stringent requirements on the control of the 
neutrino cross section. A quantitative discussion 
of this question is still not complete. For exam- 
ple, possible effects for the Minos experiment have 
been discussed by Adam Para Q . In experiments 
that measure the shape of the survival probabil- 
ity -P( 1 y tI ^i/ M )(-E'i/) as a function of the neutrino 
energy E v (for a fixed pathlength L) it is essen- 
tial to have a good reconstruction of the neutrino 
energy. This also depend to a good understand- 
ing of the neutrino cross section, to correct the 
visible energy for the energy absorbed in the ex- 
citation of the target nucleus, and the average 
composition and multiplicity of the final state. 
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An extremely good knowledge of the neutrino 
cross section will be required in case the planned 
neutrino factories will be constructed and used for 
the determination of CP violating effects in neu- 
trino flavor transitions. In this case one wants 
to compare the rates of processes such as: v e — > 
— * A* - to and T> e — > — > fi + and one will 
need a control of the ratio ov/ctu to a better pre- 
cision than the size of the expected asymmetries. 
An interesting contribution was given by Yoshihia 
Obayashi ||, for the JHF collaboration. The JHF 
project |t]] plans to use the 50 GeV proton beam 
of the Jaeri accelerator to produce an intense neu- 
trino beam directed at the Super-Kamiokande 
detector (with a 295 km baseline). The first phase 
is directed at the determination of #13 studying 
transitions — > v e , while a second phase (with 
a more massive 10 Mt Hyper-Kamiokande detec- 
tor) is aimed at a study of CP violating effects. 
A CP violation asymmetry can be defined aspl 



.4, 



CP 



-P(tV — >■ Ve) - P{VjL — > Ve) 
Pfa -> V e )+P(v^ -^Ve) 



(10) 



The oscillation probability can be schematically 
determined as: 



P{V„ -» V e ) 



expected 



(ii) 



where N e is the observed rate of e-like events, 
e e is the detection efficiency and $cxpcctcd is the 
expected v^ fluence at the detector in the ab- 
sence of oscillations. The probability P(y^ — > V e ) 
can be obtained with a neutrino beam of oppo- 
site sign. It is clear that for the determination of 
the asymmetry parameter only the ratio cr„ e / a-p c 
is relevant. In fact, considering that also the 
determination of the expected flux depends on 
knowledge of the cross section, the asymmetry 
can be rewritten as a function of the quantity 
r a = (ov e <t v ) / (oi7 c 037 ) . A part of the system- 
atic uncertainties cancel in this ratio, however it 
is not clear if, even considering this cancellation, 
the present uncertainties are sufficiently small to 
allow the measurement without a special program 
of cross section measurements. 



"In this qualitative discussion we will neglect matter ef- 
fects. 



6. Montecarlo Codes 

In all present and future neutrino experiments 
the analysis of the data will require detailed mon- 
tecarlo calculations, therefore the problem of de- 
scribing the neutrino interactions can be divided 
into two parts of similar importance: 

• The development of a theoretical framework 
to describe the interactions. 

• The development of montecarlo algorithms 
to implement this understanding. 

Of course, since a calculation of the neutrino cross 
sections with hadrons is not possible from first 
principles, the first and most important task is 
to collect experimentally the necessary informa- 
tion. The development of montecarlo algorithms 
in a non trivial problem, that has probably been 
given too little attention until this moment. All 
the neutrino experiments that have recently col- 
lected data have used different montecarlo codes, 
written by members of the collaboration, and the 
results of the different codes have not been object 
of a sufficiently careful and detailed comparison. 
One of the successes of the NUINT01 workshop, 
has been in fact the discussion between differ- 
ent groups and the comparison between different 
codes. 

The Super-Kamiokande collaboration has used 
two different montecarlo codes, the NEUT code 
Q|, originally developed for Kamiokande, and 
the NUANCE code Q (with a progenitor devel- 
oped for 1MB). The Soudan-2 detector has used 
the code NEUGEN Q that is now developed to 
be used by MINOS. The Icarus collaboration has 
performed simulations using a code developed for 
(and tuned with the data of) the Nomad detec- 
tor at CERN H|, and a second one discussed 
in |52|| . The Opera detector has performed its 
simulations with the code JETTA developed for 
CHORUS, that gives special attention to the v T 
cross section |3l| , an additional MC for Opera was 
presented at the workshop ]5q| . The montecarlo 
code of Miniboone was discussed in |55| . This is 
not the place for a critical discussion of the re- 
sults of this comparison, and only few comments 
are possible. No major discrepancy between the 
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different codes has been found, however several 
non-trivial effects have been found, and in fact 
even non negligible bugs have already been found 
and corrected for. Most codes tend to use the 
same theoretical input namely: 

(i) Llewellyn Smith |^0| for the expression of the 
quasi-elastic cross section on free nucleons. 

(ii) Smith and Moniz for the implementation of 
the Fermi gas model for quasi-elastic scattering. 

(iii) The resonance cross section of Rein and Se- 
hgal. 

(iv) The "standard" DIS formula for high W and 
Q 2 , with one of the publically available sets of 
PDF's. 

Even if the theoretical input is the same, there are 
non-trivial differences that are present, in partic- 
ular about the implementation of the Fermi-gas 
model, the "joining" of the resonance production 
and DIS scattering regimes, and the hadroniza- 
tion of the final hadronic state. The nuclear 
rescattering effects are treated (when they are) in 
significantly different ways, so effects such as pion 
reabsorption, pion elastic and charge-exchange 
scattering are very different. The field of neu- 
trino physics would certainly benefit significantly 
if some "standard" codes implementing correctly 
algorithms recognized as a possible good approxi- 
mations for the problem, could become publically 
available, as it is the case for well established 
event generators in other sub-fields of particle 
physics. The systematic uncertainties in the neu- 
trino interaction properties are still sufficiently 
large that it is probably neither possible nor re- 
ally desirable to have a unique "universal code" , 
but in the era of precision studies in neutrino 
physics, when the requirements on the event gen- 
erators will be more stringent, the development of 
more sophisticated and better tested instruments 
is clearly required. 

7. Outlook 

The study of atmospheric neutrinos with en- 
ergies around 1 GeV has yielded remarkable re- 
sults on the neutrino mass matrix. Uncertainties 
on the the neutrino-nucleus cross section in this 
energy range are of order ~ 20%. These uncer- 
tainties very likely have a small impact on the in- 



terpretation of the data, however for future stud- 
ies using accelerator neutrinos in approximately 
the same range of energy, a more precise knowl- 
edge of the cross section would be valuable, and in 
fact it will be required for studies that attempt to 
detect and measure CP violation effects in neu- 
trino oscillations. The size of the uncertainties 
on the v interaction properties can be very likely 
reduced incorporating in our description of the 
neutrino cross sections the understanding of nu- 
clear effects, hadronic resonance production, and 
the behaviour of the nucleon structure functions 
at low Q 2 that has been obtained in more recent 
studies in electron scattering. 

The near detector of the K2K experiment has 
collected the largest existing sample of v inter- 
actions in the E v ~ 1 GeV region. These data 
have a great potential to clarify several impor- 
tant questions on the neutrino interaction proper- 
ties. Hopefully a sufficient precise determination 
of the normalization and energy spectrum of the 
i/-beam will allow the measurement of absolute 
values for the cross sections; in any case impor- 
tant information about the ratios of cross sections 
can be obtained. The measurement of the rate for 
the neutral process v + A — > v + A + ix° relative 
(for example) to quasi-elastic scattering, is very 
important to put limits on oscillations of standard 
neutrinos into sterile states. 

Some natural goals for future studies could be: 

• A more quantitative estimate of the impor- 
tance of uncertainties on er„ as limiting fac- 
tor for the sensitivity of future experiments 
on neutrino oscillations. 

• A complete analysis of the K2K near detec- 
tor data to obtain the most of these remark- 
able results. 

• The development of a treatment of nu- 
clear effects that incorporates modern un- 
derstanding and calculation methods, going 
beyond the Fermi gas approximation. 

• More accurate modeling of pion and nucleon 
rescattering in nuclear matter. 

• An improved treatment of hadronic reso- 
nance production, with a smooth joining of 
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this kinematical region to the deep inelas- 
tic region, using the insight allowed by the 
electron-scattering results. 

• Development and testing of the montecarlo 
algorithms that implement our knowledge 
of the v interaction properties. 

• Last but not least: the proposal of new 
experimental studies aimed at the mea- 
surement of the most important /interesting 
properties of neutrino interactions, using 
the new intense z/-beams. See for example 
the contribution of Jorge Morfin in these 
Proceedings fl6l| . 

Neutrinos can be a precious window on the 
physics of the unification (beyond the standard 
model), but on the other hand they can also be 
an extraordinarily valuable probe on the prop- 
erties of QCD and the structure of the nucleons. 
The (V — A) structure of the neutrino interactions 
allows to explore with an axial probe the proper- 
ties of a bound system of quarks in ways that are 
very difficult to obtain in e-scattering. The cal- 
culation of the nucleon structure functions and 
their evolution in the low Q 2 region remains an 
important challenge for future QCD studies. In 
summary the measurement of the v and V interac- 
tion properties can be seen both as an instrument 
necessary to determine with high precision the v 
masses and mixing and as a a important scientific 
goal per se. In the golden era of neutrino physics, 
the motivations are powerful, the difficulties chal- 
lenging and the opportunities rich e varied. 
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